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Abstract: A compact isotope ratio laser spectrometry (IRLS) instrument was developed 
for simultaneous measurements of the D/H, 18 0/ 16 0 and 17 0/ 16 0 isotope ratios in water by 
laser absorption spectroscopy at 2.73 um. Special attention is paid to the spectral data 
processing and implementation of a Kalman adaptive filtering to improve the measurement 
precision. Reduction of up to 3-fold in standard deviation in isotope ratio determination 
was obtained by the use of a Fourier filtering to remove undulation structure from spectrum 
baseline. Application of Kalman filtering enables isotope ratio measurement at 1 s time 
intervals with a precision (<l%o) better than that obtained by conventional 30 s averaging, 
while maintaining a fast system response. The implementation of the filter is described in 
detail and its effects on the accuracy and the precision of the isotope ratio measurements 
are investigated. 
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1. Introduction 

The stable isotopes of water (in the vapor or liquid phase) are powerful tracers for the study of the 
hydrological cycle, climate change, ecological process and paleoclimatic archives (see, e.g., [1-4]), 
as well as for biomedicine [5]. Analysis of the stable isotope abundance has been the domain of 
stable isotope ratio mass spectrometry (IRMS), which generally achieves a very high precision, but 
IRMS is incompatible with condensable gases like water. In the case of water isotopologue analysis, 
reduction of water to H2 and equilibration with CO2 is commonly used for the determination of the 
2 H (or D) and 18 0 isotope ratios [6,7]. In addition, as the 17 0 12 C 16 0 and 16 0 13 C 16 0 molecules appear 
in the same mass channel, accurate and direct determination of the isotope ratio H2 17 0/H2 16 0 is 

17 

practically impossible. This ratio is usually inferred from the natural relation between the O and 

1 8 

O abundance ratios [8] or determined on O2 after fluorination of the water sample [9]. No reliable 
high precision mass spectrometric method is available for direct isotopic analysis of water. 
Moreover, since the equipment is quite bulky it cannot be used in the field and is generally incapable 
of performing real-time measurements. 

Measurements of the stable isotope ratio by optical spectroscopy, in particular by laser absorption 
spectroscopy (LAS), have attracted a growing interest in recent years [10,11]. Laser spectroscopy 
provides an excellent opportunity to perform in situ real-time continuous measurements without 
requiring chemical conversion. Nowadays isotope ratio laser spectrometry (IRLS) becomes a serious 
competitor to IRMS. Water isotope spectrometers operating in the near infrared around 1.4 urn [12-14], 
in the mid infrared at 2.7 um [15,16] and 6.7 um [17,18] have been developed. Cavity enhanced 
spectroscopy approaches have been implemented for achieving high sensitivity and high precision [18-20]. 
High precision IRLS instruments are now commercially available from LGR Inc. (Mountain View, CA, 
USA) (using off-axis Integrated Cavity Output Spectroscopy, known as OA-ICOS) [21-23] and from 
Picarro Inc. (Santa Clara, CA, USA) (using Cavity Ring Down Spectroscopy— CRDS) [24-27]. 

In this paper, we report on the development of a compact IRLS instrument for simultaneous real-time 
measurements of the D/H, 18 0/ 16 0 and 17 0/ 16 0 isotope ratios in water by laser absorption spectroscopy 
at 2.73 um. The 2.7 um fundamental stretching mode is about one order of magnitude stronger than the 
band at 1.4 um, and even the bending mode near 6.7 um is not much stronger than that the 2.7 um 
band. More important is that the selected absorption lines at 2.73 jam for measurements of water 
isotopologues (H2 18 0, H2 16 0, H2 17 0 and HDO) may have relatively similar line strengths and similar 
ground state energies (Table 1). Recently available, commercial 2.73 um distributed feed-back (DFB) 
lasers offer the following advantages compared to a color center laser previously used for water 
isotopologue measurement [15]: single mode and single frequency emission at room temperature, high 
electronic bandwidth (as high as kilohertz repetition rates), in combination with compactness (mounted 
inside a TO-8 package) and cost-effectiveness [16]. 

In the present work, special attention is paid to the spectral data processing by use of digital 
filtering techniques to improve the measurement precision. Oscillatory structure on spectrum baseline 
affecting the measurement precision was analyzed. Fourier analysis of spectral residuals resulting from 
a fit was performed, which enabled the identification and filtering of the noise in spectral signals with 
the help of Fourier filter. The paper also provides a detailed description of Kalman filtering of the data, 
as recently introduced by us to the field of isotope ratio measurements [28]. We demonstrate that a faster 
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temporal response (Is) can be achieved by Kalman filtering with a precision better than that obtained 
by conventional 30 s averaging. Effects of Kalman filtering on the accuracy and the precision of the 
isotope ratio measurements are also studied in the present work. 

2. Experimental Consideration 

2.1. Isotopic Composition Analysis by Laser Absorption Spectroscopy 

According to the Beer-Lambert law of linear absorption, the absorbance A(v) can be related to the 
transmitted light I(v) and incident intensity I 0 (v) intensities: 

A(v) = In (/ 0 (v)//(v)) = Ca{v)L (1) 

where C is the number density of absorbing species (in mol/cm ), a(v) is the frequency-dependent 
absorption cross section (in cm /mol), and L is the optical absorption path length (in cm). The 
integrated absorbance A/ (in cm -1 ) can be written as: 

Aj = j A(v)dv = j ln(/ 0 (v)) / I{v))dv = CL\ a(v)dv = CLS(T) I n (2) 

with 5(7) the temperature-dependent molecular line absorption intensity in cm _1 /(mol-cm~ 2 ) and n the 
fractional isotope abundance; both can be found in the HITRAN database [29]. The dependence of the 
absorption line intensity 5(7) on the temperature T can be described as follows [13]: 



S(T) = S(T Q ) 



f T ^3/2 

1 0 



exp 



hcE 0 



(3) 



with To = 296 K, Eo the lower level energy expressed in cm -1 , h the Planck constant (J s), c the speed 
of light in vacuum (cm/s) and k the Boltzmann constant (J/K). 

The isotope ratio can thus be determined from the ratio of the integrated areas A/ and the absorption 
line intensities S of the major and minor isotopic components: 

ax ca , a 

R = [C x )/[C a ) = ^-x^- (4) 
a; S x ln x K ) 

where x refers to the rare isotopic species (H 2 O, H 2 O, or HDO), a represents the abundant isotopic 
component (H2 16 0), and R represents the ratio of the rare to the abundant isotopologues abundances. 

The relative deviation of the isotope ratio in water with respect to the international standard 
reference known as Vienna Standard Mean Ocean Water (VSMOW), is expressed in terms of the S-value: 

^/ oo ),^™^xl000 (5) 

RvSMOW 

1 8 

where, by international convention, R vsm0 w takes on the values: R vsm0 w = 0.0020052 for O, 0.0003799 

17 2 

for O and 0.00015576 for H [30]. It is noted, however, that substitution of Equation (4) into 
Equation (5) results in an expression for the delta-value in which the ratios of the line strengths cancel, 
provided that the sample and reference spectra are recorded at the same temperature (see, e.g., 
Reference [10]). In practice, the laser spectrometer is calibrated against working reference standards 
with isotopic ratios that are well known on the international two-point scale (determined by the 
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consensus values for VSMOW and Standard Light Antarctic Precipitation — SLAP), recommended by 
the IAEA [30]. A few years ago, the International Atomic Energy Agency (IAEA) produced the new 
reference material VSMOW2 to replace the exhausted reference material VSMOW. Repeated 
measurements of VSMOW and VSMOW2 by a number of international laboratories have shown that 

1817 2 1817 

they have indistinguishable 8 O, 8 O and 8 H values, and thus indistinguishable R vsmow for O, O 
and 2 H [31]. 

In the present work, the laser instrument was calibrated against the working standards GS-49 
(8 18 0 = 0.39%o, 8 17 0 = 0.21% o , and S 2 H = 1.7%o, with respect to VSMOW, determined by repeated 
IRMS analyses at the Center for Isotope Research of the University of Groningen). The accuracy of the 
laser instrument was evaluated by measurement of another working standard GS-42 with different 
isotopic composition (8 18 0 = -24.62%„, 8 17 0 = -13.1%o, and 8 2 H = -187 .7%o). Bottled water 
(Vittel, France) was used as unknown sample material. 

2.2. Selection of the Absorption Lines 

Selection of suitable absorption lines for water isotopic ratio measurements is one of the most 
important aspects in instrumental design since the choice of absorption lines has a direct impact on the 
instrumental performance in terms of measurement sensitivity, precision and selectivity. 



Table 1. Comparisons of line intensities, ground state energies and temperature 
coefficients between different IRLS operating at different wavelengths in the infrared 
spectral region. 



Reference 


Isotopologue 


Frequency 

(cm" 1 ) 


Intensity 
(10 23 cm-mol 1 ) 


Ground state 
energy (cm -1 ) 


Temp. coef. 
at 296 K (K 1 ) 




H 18 OH 


7183.5858 


0.62 


733.7 


6.9%o 


[10,11] 


H 16 OH 


7183.6858 


0.31 


661.6 


5.7%o 


H 17 OH 


7183.7354 


0.12 


95 


-3.4%o 




H 16 OD 


7183.9727 


0.035 


156.4 


-2.5%o 




H 18 OH 


1483.9261 


8.4 


550.5 


4.1% 0 




H 16 OD 


1484.1064 


2.3 


225.9 


-1.3°/oo 


[16] 


H 16 OH 


1484.2573 


1.8 


1899.0 


26.2%o 


H 17 OH 


1484.5109 


2.0 


205.5 


-1.6°/oo 




H 18 OH 


1484.9716 


10.0 


325.2 


0.3%o 




H 16 OH 


1485.1336 


6.2 


1907.6 


26.4%o 




H 18 OH 


3662.9196 


2.1 


398.3 


1.5% 0 


[13,14] and 


H 16 OH 


3663.0452 


8.5 


586.4 


4.6%o 


our work 


H 17 OH 


3663.3213 


7.2 


224.3 


-1.4°/oo 




H 16 OD 


3663.8419 


1.2 


100.4 


-3.4%o 



18 17 2 

Measurements of the isotopic ratios 8 O, 8 O, 8 H of the stable isotopologues of water require 
probing of absorption lines of the four isotopologues H 2 O, H 2 O, H 2 O and HDO within the laser 
tuning range. It is very desirable that the used absorption lines exhibit similar absorption depths at natural 
abundance with an absorption intensity as large as possible (for high sensitivity and high precision 
measurements), are free from interference from the same or other species (high selectivity consideration), 
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and have similar ground state energy (in order to minimize the effects of temperature-dependent line 
intensities). The spectral region near 2.73 um, covered by recently commercially available DFB lasers, 
meets these requirements to a high degree. Parameters of the four molecular lines selected in the 
present work are summarized in Table 1. Parameters of the previously used lines at 1.4 um and 6.7 um 
are also given for comparison. 

The line absorption strength depends on the population of the ground state level and this population in 
turn depends on the temperature. Temperature drift during measurements may introduce a systematic 
error in the isotope ratio determination. The temperature coefficients defined as AS/S(To) for the 
selected lines are listed in Table 1 . A thermal drift of 1 K would lead to a relative variation in the line 
strength of +1.5%o, +4.6%o, -1.4%„, -3.4%o for the selected lines of H 18 OH, H 16 OH, H 17 OH and 
H 16 OD respectively. Therefore, realization of a temperature stability better than 0.1 K is essential for 
high precision measurements. 

3. Experimental description 

3.1. Experimental Set-Up 

The experimental set-up, mounted on a 50 x 70 cm optical breadboard, is depicted in Figure 1. The 
laser source was a room temperature single mode DFB diode laser operating at 2.73 um (Nanoplus 
GmbH, Gerbrunn, Germany). It is tunable from 2729 to 2732 nm with an output power of 2 mW. 
A Thorlabs (Newton, NJ, USA) ITC 502 diode laser controller provided the laser temperature control 
and laser drive current. The diverging laser beam was collected by an off-axis parabolic mirror PMl with 
an effective focal length (EFL) of 25 mm. The laser beam was then transformed into a quasi-parallel beam 
with a diameter of ~4 mm by a combination of an antireflection coated CaF2 lens Fl (f = 200 mm) and 
an off-axis parabolic mirror PM2 (EFL = 50 mm). 

Figure 1. Optical layout. PMl and PM2: parabolic mirrors with effective focal length of 
25 mm and 50 mm, respectively; Fl and F2: lenses of focal length of 200 mm and 50 mm, 
respectively. 



Multi-pass cell Beam shutter 




F2 



The collimated beam was then divided into two parts. The first part (-8%) was reflected by a beam 
splitter (CaF2) and directed to a homemade Fabry-Perot etalon for spectral metrology. The frequency scale 
was linearized by means of the interference fringes produced by the etalon consisting of two air-separated 
uncoated CaF 2 plates with a free spectral range (FSR) of -0.0283 cm -1 . Positions of the H 2 0 vapor 
absorption lines provided by the HITRAN database [29] were used as absolute frequency reference for 
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laser frequency calibration. The main part of the laser beam was coupled into a Herriott cell with an 
absorption path of 20 m in a 36-pass configuration. The emerging absorption signal was focused with 
a 50 mm CaF 2 lens F2 onto a LN 2 cooled HgCdTe detector (J15D22-M204-S01M-60). A home-built 
bridge circuit was employed to realize DC coupling of the HgCdTe detector to a low noise 
preamplifier (Model 5113, EG&G, Albuquerque, NM, USA). 

The signal background (corresponding to "laser off), determined by the dark current level of the 
detector, was found to be fluctuating by about 1% over a 1 h time interval. This probably was due to 
variation of the detector temperature or its DC power supply. In order to correctly retrieve absorption 
spectra, a beam shutter (Thorlabs, SH05) was placed before the detector for background level 
acquisition at the beginning of each absorption spectral scan. 

The laser frequency was periodically scanned at a rate of 10 Hz across the absorption lines of the 
water isotopologues H 16 OH, H 17 OH, H 18 OH, and HDO by means of a triangular wave voltage. The 
output signal from the preamplifier was digitized with a laptop using a 16-bit analogue/digital data 
acquisition card (DAQ Card-6036E, National Instruments, (Austin, TX, USA) controlled with a Lab 
Windows-based program. 

3.2. Experimental Protocol 

Careful attention has been paid in the present work to sample mass effects and sample memory 
effects that affect the determination of the isotopic ratios, as discussed in detail by, among others, 
Kerstel et al. [15] and more recently Lis et al. [21]. Liquid water samples of 12 uL were injected into 
the pre-evacuated gas cell through a silicon membrane using a syringe resulting in a water vapor 
pressure of -4.5 mbar inside the gas cell (the saturated vapor pressure is -42 mbar), which corresponds 

17 3 

to a H 2 0 molecule number density of 1.1 x 10 mol/cm in the cell. Besides temperature fluctuations, 
variations in sample pressure inside the cell, resulting from sample injection by a syringe through 
a silicon membrane, will have a strong effect on the precision and accuracy of the measurements. 
Although constant volume was used for both sample and reference standard injections in our 
experiments, variances in the H 2 0 number density in the cell may occur due to H 2 0 leakage through 
the silicon membrane, over-tightening or under-tightening of the silicon membrane. In order to 
minimize the impact of this effect on the final results, only those injections resulting in a gas cell 
pressure within +0.1 mbar of 4.5 mbar were accepted. Furthermore, in order to avoid sample memory 
effects due to the "stickiness" of water on the gas cell wall, the first three injections were discarded for 
each isotope ratio determination. 

The temperature of the gas cell was actively controlled to 30 °C by the use of a heater band, and 
maintained constant within +0.1 °C by a PID controller. The cell temperature was monitored with 
calibrated platinum resistors (PtlOO) with an accuracy of 0.03 °C and a precision of 0.01 °C. 
No temperature gradient along the cell axis was observed within the measurement precision of the 
temperature sensors. 

In our experiment, a slight drift of the laser wavelength with time has been observed. In order to 
minimize the effects of the instrumental drift on the determination of A/, only 10 laser scans were 
co- added, resulting in a raw 8- value at 1Hz averaged data acquisition rate. Further improvement in 
precision has been achieved through the use of Kalman filtering technique. As can be seen in the 
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results presented in Section 5, the use of Kalman filtering approach allows minimizing the effect of 
laser wavelength drift since it filters the measurement data at 1 s time intervals. 

4. Data Processing and Retrievals 

Figure 2 shows a typical experimentally recorded spectrum of H 2 0 isotopologues around 2.73 um. 
The spectrum, resulting from an average of 10 laser scans, was recorded at a pressure of 4.5 mbar and 
a temperature of 30 °C. Spectral data processing for isotopic composition determination is discussed in 
detail in the following subsections. 

Figure 2. Experimental spectrum of H2O isotopologues around 2.73 um. At the bottom is 
shown the fringe from an etalon that was used for frequency metrology in combination 
with the H 2 0 line positions provided by the HITRAN database. 
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4.1. Spectral Fitting Algorithm 

In order to determine the integrated absorbance A/ (the area under the absorption line profile), 
absorption spectra were fitted to Voigt [32] and Galatry [33] profiles, using a Levenberg-Marquardt 
multi-line fitting algorithm. In the fit procedure, the baseline approximated with a 4th-order 
polynomial was experimentally found to well represent laser power ramp, allowing good removal of 
the variation of laser power in the fitted residual. Apart from the four selected lines of H 18 OH, H 16 OH, 
H 17 OH and H 16 OD, two H 16 OD lines at 3663.2250 cm" 1 and 3663.3879 cm" 1 on either side of the 

17 

Ff 'OH line were also taken into account in the spectral fitting program. As indicated by the fit residuals 
shown in Figure 3, the Galatry profile (b) fit resulted in a better residual than using the Voigt profile (a). 
In fact, the Voigt profile does not take into account correlations between molecular velocities and 
collisional processes. Corrections for velocity-changing collisions are included in the soft-collision 
(Galatry) model [33]. Speed dependence of the relaxation rates can be accounted for with the 
speed-dependent (SD) Voigt profile [34]. Both Galatry and SD- Voigt profiles are reported for 
recovering the experimentally observed lineshapes, and Galatry profile essentially behaves like the 
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speed-dependent Voigt model [35]. However, Galatry profile and SD-Voigt profile are based on 
completely different processes. As can be seen in Figure 3c, d, relatively big residuals around absorption 
peaks are still evidenced even using Galatry model, which leads to a bigger uncertainty in the 
determination of the integrated area Ai. In order to significantly minimize the uncertainty in the 
determination of the integrated area under the absorption features, a more sophisticated model could 
be adopted in the future (for instance, the speed-dependent Galatry profile), to take into account the 
narrowing due to the speed-dependence of relaxation rates as well as the averaging effect of 
velocity-changing collisions. 



Figure 3. H2O isotopologue absorption spectrum measured in this work (raw data). 
Spectral data were fitted to Voigt and Galatry profiles. Residuals resulting from each fit are 
shown: (a) from Voigt fit and (b) from Galatry fit, respectively, (c) shows the fit results 
using Galatry profiles in combination with a Fourier series function (FS) to remove 
oscillation structure from the spectrum baseline, (d) shows a fit residual from a spectrum in 
the absence of an undulating structure on the baseline. 
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4.2. Suppression of Oscillation Structure on the Baseline 



A periodic oscillatory structure on the baseline was observed in both the Voigt (Figure 3a) and the 
Galatry fit residuals (Figure 3b). This kind of sinewave-like undulation, usually assumed to be fringes 
resulting from optical interference, was, however, also observed with laser turned off. This means that 
the periodic oscillatory structure was caused by an electrical perturbation. In order to further improve 
the measurement precision, we applied a Fourier filtering approach [18] to suppress this oscillatory 
baseline structure. For this purpose, Fourier transformation of the residuals of the fit was performed in 
order to determine the oscillation frequencies. As can be seen in the noise power spectral density graph 
(Figure 4), the noise at frequencies higher than 0.001 Hz was associated with the Voigt fit residual as it 
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has been almost completely removed in the Galatry fit. We determined that the noise peaks at 
frequencies lower than 0.001 Hz were associated with the baseline oscillatory structure. In order to 
take this undulation structure into account, the baseline was then modeled as a composition of a 4th-order 
polynomial and a Fourier series (FS) function of sine and cosine waves (for simulation of the 
sinewave-like oscillatory structure of the baseline). The transmitted power signal recorded by 
experiment was thus fitted to the following function: 



where m is acquired data point number, Po(m) is a 4th-order polynomial representing the laser power 
ramp, Ak and Bk are the oscillation amplitudes, and a>k is the oscillation frequency determined from the 
Fourier transform. In our experiments, only the most prominent three frequencies have been taken into 
account: coi = 3.662 x 10" 4 Hz, co 2 = 4.882 x 10" 4 Hz and co 3 = 7.324 x 10" 4 Hz. Gi(m) is the Galatry 
function which describes the direct absorption line shape for the H 18 OH (3662.9196 cm -1 ), H 16 OH 
(3663.0452 cm" 1 ), H 16 OD (3663.2250 cm" 1 ), H 17 OH (3663.3213 cm" 1 ), H 16 OD (3663.3879 cm" 1 ) and 
H 16 OD (3663.8419 cm" 1 ) lines. The Galatry profile can be defined by: 



frequency, y D (cm ! ) the Doppler half-width, y L (cm ) the collisional half-width, and /? c (cm l ) 
representing the average effect of collisions on Doppler broadening. 

Figure 4. Discrete Fourier transform of the fit residuals resulting from different spectral 
fitting algorithms. 



P = P 0 (m) exp (-G(m))(l + £ (A k sm(co k ) + B k cos{w k ))) 



(6) 




(7) 
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To obtain the integrated absorbance A/, we fitted the F^ function (Equation (6)) to the 
experimentally recorded data with fixed isotope-dependent Doppler width y D (scaled with the square 
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root of the isotopologue mass), while the following parameters were determined in the fit: the 
integrated absorbance Aj, the line center position vo, the Lorentzian pressure broadening width y L , the 
Dicke narrowing parameter /? c , and the baseline modeling parameters. 

Figure 3c and the black lines in Figure 4 show the results after suppression of the periodic undulation 
structure by fitting the spectrum to Equation (6). As can be observed in Figure 5, compared to the results 
obtained with a simple 4th-order polynomial as baseline model (left panel in Figure 5), baseline correction 
using a composition of a 4th-order polynomial and a Fourier series function (central panel in Figure 5) 
leads to a reduction of up to 3-fold in the standard deviation (SD) in the isotope ratio determination. 

Figure 5. The 18 0/ 16 0, 17 0/ 16 0 and D/H isotope ratios obtained with different baseline 
modelings. Left panel: baseline described with a 4th-order polynomial. Central panel: 
baseline modeled using a 4th-order polynomial and a Fourier series function (with three main 
frequency components). Right panel: no oscillation structure in original spectrum baseline. 
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Though application of Fourier filtering can improve the measurement precision by removing 
oscillatory structure from the baseline, it is hard to completely account for the exact baseline structure 
with only three Fourier components. In fact, the residuals of Figure 3c still show signs of a residual 
periodic structure. We therefore proceeded to eliminate the source of the electrical perturbation, which 
was determined to be related to an electronic component failure in the preamplifier. In the case of no 
electrical perturbation, as shown in Figure 3d, a 3.7-fold improvement in the measurement precision 
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was obtained (right panel in Figure 5), in comparison with the results achieved using a Fourier filtering 
(middle panel in Figure 5). 

4.3. Allan Variance 

The measurement precision is usually affected by the instrument instability and measurement errors 
related to sample handing and injection (e.g., incomplete evacuation of the gas cell between two 
consecutive measurements will lead to memory effects and thus affect the measurement precision). 
Another important limiting factor to achieving high precision is the signal-to-noise ratio (SNR) of the 
spectral data. In LAS, high SNR can be obtained by enhancing the signal (by selecting stronger 
absorption lines and using long absorption path length or cavity enhanced spectroscopy) and reducing 
the noise (by averaging N laser scans and using modulation or/and balanced-beam detection 
techniques). With the signal averaging approach, the optimal averaging number N, limited by the 
stability of the instrument, can be determined by an Allan variance analysis [36]. Alternatively, the 
individual spectral scans can be processed, and the resulting isotope ratios may be further averaged 
(within the system's stability time) to obtain the desired precision level. In the present work, due to 
laser wavelength drift, the maximum averaging number of the laser scans was limited at about N = 10, 
corresponding to an acquisition time of 1 s. The 1 s raw data were processed to provide 1 s raw 8 values 
with a la precision of 7.8%o for 5 18 0, 6.6%o for 5 17 0, and 8.0%o for 8 2 H of a bottled water (Vittel, France). 

Figure 6. Measurement results of a bottled water (Vittel, France). The upper three panels 
show raw measurement of the S-value for 18 0, 17 0, 2 H (dots) and the corresponding 
Kalman filter output for a q- value of 800 (lines). The Allan variance plotted in the lower 
panel shows an optimal averaging time of -30 s for the present IRLS system. 
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In order to further improve the measurement precision, the 1 s raw 8 values are further averaged. 
Allan variance analysis has been performed to determine the optimum averaging number. As can be 
seen in Figure 6 (lower panel), the optimum averaging time for the present instrument was -30 s. We 
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then used conventional averaging of 30 measured S-values corresponding to the optimum averaging 

1 Q IT 

time of 30 s. Herewith the precision was improved to 1.4%o for 8 O, l.l%o for 8 O, and 1.5%o for 
8 H, respectively. 

4.4. Kalman Filtering 

Though signal averaging enables one to improve the measurement precision [1,11], such 
"post-processing" results in a slow temporal response of the system. For certain specific applications, such 
as the on-line monitoring of exhaled breath [37], it is highly desirable to be able to perform real-time 
measurements with high sensitivity and precision, while maintaining a fast system response. Kalman 
filtering, being an adaptive filtering technique uses a recursive procedure for "true value" prediction 
based on the previously determined values. It can efficiently remove the shot-to-shot variability related 
to the real-time noise in the measured data with minimal deformation of the physical quantity to be 
measured. Kalman filtering has been successfully applied before to real-time trace gas concentration 
measurements [38,39]. We recently introduced this technique to the field of isotope ratio 
measurements [28] in order to perform fast and high precision measurements. In this section, we 
describe in detail the theoretical consideration of a Kalman filter model for application to isotope 
ratio measurements. 

Using a linear stochastic difference model, the true isotope ratio S k+1 at time k+l is evolved from 
the value 8 k given at k according to: 

4+i=4 +w * (8) 

At time k the measured isotope ratio Zk of the true value S k can be expressed as follows: 

Z k =S k +v k (9) 
where Wk and Vk are uncorrected random variables related to the process variability and the 

2 2 

measurement noise with corresponding covariance of c w and a v respectively. Note that a hat above 
a variable indicates an estimated or predicted quantity, and a superscript negative sign above a variable 
indicates an a priori quantity. The recursive procedure involved in Kalman filtering may be considered 
to consist of two parts: the "time update" and the "measurement update". 

The time update procedure projects forward in time current isotope ratio and error variance estimates 
to obtain a priori estimates for the next time step. The prediction equations can be expressed as: 

4+1=4 (10) 

Pk + i=P k +^l (11) 

The measurement update incorporates a new measurement into the a priori isotope ratio estimate to 
obtain an improved a posteriori estimate as filtered isotope ratio value, 8. The filtering equations 
involved in the measurement update are: 

K k =p-{P- + C7 2 v r (12) 
4=4 + ^(^-4) (13) 
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P k =P k (X~K k ) (14) 

where 8~ k+1 is the predicted isotope ratio estimate at step k+1, 8 k is the filtered isotope ratio value 
estimated at step k. P~ k+ i is a priori estimate error variance, P k is a posteriori estimate error variance, 
K k is the Kalman gain that weights the measurement residual, defined as the difference between an actual 
isotope ratio measurement Zk and a measurement prediction S~ k as shown in Equation (13). 

In practical applications, the measurement noise a v and the true § variability a 2 M , (related to real 
isotope abundance variation and real-time drifts resulting from laser frequency shift, thermal 
fluctuation, pressure variation, etc.) should be well defined. 

Information on the measurement noise variance o 2 v is usually known, because it depends on the 
quality of the measurement instrument, while the variance c w of the true isotopic ratio variability is 

2 2 

quite subjective. Whereas both a v and a » vary with real variations of the isotope abundance, once 

2 2 

the measurement system has reached equilibrium, the ratio of c v to a w should be constant. We 
therefore define: 

q = o 2 v /o 2 w (15) 



Figure 7. (a) Plots of the standard deviation in water (Vittel) isotope ratio determination 
in function of the q values involved in Kalman filtering; (b) The relative deviation in 
water (GS-42) isotope ratio measurement accuracy, {8-8)1 8 , plotted in function of the 
q-parameter. 
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The parameter q is known as the filter tuning parameter. The choice of q depends on the particular 
instrument and its application environment. For a larger q value, the system's response time is longer 
to follow real-time variation in isotopic composition. Conversely, the filtering is less efficient in 
removing shot-to-shot real-time noise when a smaller q value is used. Figure 7a shows the la precision 
(standard deviations) in isotope ratios determination of a bottled water (Vittel) as a function of q. As 
shown in Figure 7a, the improvement in the precision is clearly evident, as the standard deviation 
decreases exponentially with increased q value. The decreasing trend for 5 18 0 is toward 
a steady state value for q > 800. For 8 17 0 and 8 2 H, the decreasing trend slowed down for q > 800, 
while toward a steady state value for q > 2000. This may be explained by small drift in 8 17 0, 8 H 
resulting a Gaussian distribution, high q value is needed to get steady state, because the larger q value, 
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the more efficient is the filtering 8-value fluctuations due to random noise. In the present work, a value 
of q = 800 was chosen as a compromise between fast temporal response and high filtering efficiency 
(thus high measurement precision) for the developed laser instrument. In the present work, a value of 
q = 800 was chosen as a compromise between fast temporal response and high filtering efficiency 
(thus high measurement precision) for the developed laser instrument. In our experiment, c v was 
determined by the 8-value variance deduced from the first 10 raw measurements, and c w was 
calculated by dividing c v by q. 

The effects of the Kalman filtering on the measurement accuracy were investigated in the present 
work. For this purpose, the working standard GS-42 with well known isotope ratios was used. The 
relative deviation in water (GS-42) isotope ratio measurement accuracy, {8-8)18 , is plotted in 

18 2 

function of the q-parameter in Figure 7b. As can be seen, the relative deviation for 8 O and 8 H are 
almost constant in the range of ±1% with increased q-value. Though the relative deviation for 8 O 
varies evidently with q-value, the values are under a range of +0.7%. The Kalman filtering is more 
effective to remove fluctuations due to random white noise. However, apart from the measurement 
random noise, temperature fluctuations of the gaseous sample and baseline drift may result in drift 
over time in 8 values, which cannot be removed by the Kalman filtering. Because of slow drift in 8 17 0 
values, the fluctuation of 8 17 0 mean for different q-value is almost two times larger than 8 18 0 and 8 2 H. 
The maximal differences in relative deviation for the measurement accuracy are 1.6 x 10" 3 for 8 18 0, 
1 x 10" 2 for 8 17 0, and 5 x 10" 4 for 8 2 H in Figure 7b. 

5. Results and Discussion 

Table 2 summarizes the measurement accuracy of our IRLS instrument, compared to the IRMS 
calibrated value of GS-42. 



Table 2. IRLS vs. IRMS: measurement accuracy comparison of isotope ratios in water 

1 /9 

reference GS-42. Where SE = SD/N 1 , is the standard error of the mean values for N 
measurements. IRLS results were obtained using 1 s Kalman filtering with q = 800. 





N 


5 ls O (SE) 


5 17 0 (SE) 


5 2 H (SE) 


IRMS (CIO) 


3 


-24.62 (0.02) 


-13.1 (0.1) 


-187.7 (0.3) 


IRLS (present work) 


5 


-24.83 (0.50) 


-11.7 (0.3) 


-186.2 (0.6) 



The deviation of the measured values from the reference was -0.2%o for 8 18 0, 1.4%o forS 17 0 and 
1.5%o for 8 2 H. The large measurement deviation, especially for 8 17 0 and 8 2 H, would be mainly caused 
sample mass effects. Though careful attention has been paid to sample amount effects in the present 
work, more precise control of injection volume and sample pressure was still needed. The la standard 
deviation is 1.2%o for 8 18 0, 0.7%„ for 8 17 0, and 1.4%o for 8 2 H for five injections of water (GS-42). 
Moreover, the precision of 0.8%o for 8 18 0, 0.6%o for 8 17 0 and 0.9%o for 8 2 H with 1 s Kalman filtering 
also had certain influence on the measurement deviation. It should be noted that the S 17 0 value in 
water could not be directly determined with IRMS; the value given in the table is not a measured 
value, but was calculated based on a natural relation between the O and O abundance ratios [8]: 
S 17 0 = (1+ S 18 0) a5281 - 1. 



Sensors 2014, 14 



9041 



The measurement precisions using 1 s Kalman filtering (q = 800) are given in parentheses in Table 2 
in terms of the standard error of the mean values, SE (= SD/N 1/2 ). We converted the standard deviation 
into SE with the same N values for the reason of comparison with the standard working reference GS-42. 

Figure 6 shows raw measured results of a bottled water Vittel (black dots) in comparison with 
Kalman filtering results (black lines). With a q-value of 800, the Kalman filtered data obtained in 1 s 
reached a precision level better than the results obtained from conventional 30 s averaging. The la 
standard deviation has been improved to 0.83%o for 8 18 0, 0.57%„ for 8 17 0, and 0.91%o for S 2 H, with 
a measurement time of 1 s. Comparison of the measurement precisions is given in Table 3.1s Kalman 
filtering 8 values (with q = 800) are compared with the raw measured 8 values from the average of 10 
laser scans (Is) and the data from conventional averaging of 30 measured 8-value within the optimal 
averaging time of 30 s. 

Table 3. Measurement precision comparison: Raw measured 8-value from the average of 
10 laser scans in 1 s, data from conventional averaging of 30 measured 8-values within the 
optimal average time (30 s) determined by an Allan variance analysis, and 1 s Kalman 
filtering data with q = 800. 



Method / Time 


5 ls O (SD) 


Measurement precision 
5 n O (SD) 


5 2 H (SD) 


Raw measurement /Is 


7.8 %o 


6.6 %o 


8.0 %o 


Averaging 30-5 / 30 s 


1.4 %o 


1.1 %0 


1.5% 


Kalman filtering /Is 


0.8 %o 


0.6 %o 


0.9 %o 



It is also worth noting, as can be seen in Figure 6, that the Kalman filtering does not affect the mean 
value of the data, and therefore does not affect the accuracy of the measurement. 

6. Conclusions 

We have described a compact IRLS instrument for simultaneous measurements of the D/H and 
18 0/ 16 0, 17 0/ 16 0 isotope ratios in water by laser absorption spectroscopy at 2.73 um. We demonstrated 
the potential of application of Kalman filtering for fast and high precision isotope ratio measurements. 
The measurement precision achieved at 1 s Kalman filtering time intervals is better than that obtained 
by conventional 30 s averaging. The Kalman filter can be optimized to filter out the maximum amount 
of shot-to-shot real-time noise while following the real variation in measured physical quantity. The 
impact of the Kalman filtering on the measurement precision and accuracy were investigated. The 
measurement precision was improved clearly with high q value, while measurement accuracy was less 
impact by the Kalman filtering and mainly determined by calibration with known standard materials. 
The Voigt and Galatry profiles were used to investigate the possible influence of the choice of the line 
shape model on the determination of the integrated absorbance. The Galatry profile fit resulted in 
a smaller residual than using the Voigt profile. However, relatively big residuals around absorption 
peaks are still evidenced even using Galatry model, which leads to a bigger uncertainty in the 
determination of the integrated area Ai. A more sophisticated model should be adopted in the future 
(for instance, the speed-dependent Galatry profile), to take into account simultaneously the narrowing 
due to the speed-dependence of relaxation rates and to the averaging effect of velocity-changing 
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collisions. Fourier analysis of fitted spectral residuals allowed for determination of the Fourier 
frequency components in the undulation structure on the baseline. Reduction of up to three times in 
standard deviation in the isotope ratio determination was obtained by application of Fourier filtering to 
remove the undulation structure from spectrum baseline. However, a further reduction by a factor of 3.7 
was obtained by the elimination of the electrical noise source. Measurement accuracy comparison of 
isotope ratios in water reference GS-42 between IRLS and IRMS was performed. Measured value of 

18 17 2 

8 O was more close to reference value, while measured values 8 O and 8 H presented a deviation 
from reference value. The measurement deviation would be mainly caused sample mass effects. 
Though careful attention has been paid to sample amount effects in the present work, more precise 

1 Q 

control of injection volume and sample pressure was still needed. Moreover, the precision of 8 O, 
8 17 0 and S 2 H with 1 s Kalman filtering also had certain influence on the measurement deviation. 

Further improvements in isotope ratio determination precision and accuracy can be envisaged as 
follows: (1) precise control of injection volume and sample pressure; (2) real-time calibration by 
alternating the introduction into the cell of a reference and the sample [21,40,41]; (3) real-time 
precise control and measurement of gas cell temperature for correction for temperature drifts during 
the measurements. 
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